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ABSTRACT 

The statistics of gravitationally lensed arcs was recognised earlier as a potentially pow- 
erful cosmological probe. However, while fully numerical models find orders of magnitude 
difference between the arc probabilities in different cosmological models, analytic models 
tend to find markedly different results. We introduce in this paper an analytic cluster lens 
model which improves upon existing analytic models in four ways. (1) We use the more re- 
alistic Navarro-Frenk- White profile instead of singular isothermal spheres, (2) we include the 
effect of cosmology on the compactness of the lenses, (3) we use elliptical instead of axi- 
ally symmetric lenses, and (4) we take the intrinsic ellipticity of sources into account. While 
these improvements to the analytic model lead to a pronounced increase of the arc probabil- 
ity, comparisons with numerical models of the same virial mass demonstrate that the analytic 
models still fall short by a substantial margin of reproducing the results obtained with nu- 
merical models. Using multipole expansions of cluster mass distributions, we show that the 
remaining discrepancy can be attributed to substructure inside clusters and tidal fields con- 
tributed by the cluster surroundings, effects that cannot reasonably and reliably be mimicked 
in analytic models. 



1 INTRODUCTION 



Many authors pointed out that the statistics of gravitationally lensed 
arcs in galaxy clusters may be a powerful tool for constraining cos- 
mological models. They are rare events caused by a highly non- 
linear effect in cluster cores, and are thus not only sensitive to the 
cosmologically highly variable number density of galaxy clusters, 
but also to their internal structure. 

The expected number of giant arcs, usually defined as 
arcs with a length-to-width ratio exceeding ten and apparent B- 
magnitude less than 22.5 (Wu & Hammer 1993), changes by orders 
of magnitude between low- and high-density universes according to 
the numerical models described in Bartelmann et al. (1998). They 
use the ray-tracing technique for studying gravitational lensing by 
galaxy cluster models taken from N-body simulations (see also 
Bartelmann & Weiss 1994; Bartelmann 1995; Bartelmann, Stein- 
metz & Weiss 1995; Meneghetti et al. 2000, 2001). This allows 
the most realistic description of the cluster lenses because all ef- 
fects which could play an important role for the lensing phenom- 
ena are by construction taken into account. Cluster asymmetries, 
substructure, and the tidal field of the surrounding matter distribu- 
tion are known to have substantial effects on arc statistics, and they 
are automatically included if the cluster simulations are suitably 
designed. 

In order to extract useful cosmological constraints from arc 
statistics, it is essential to perform simulations on a fine grid in the 
cosmological parameter space. However, given the long computa- 
tion times required for full numerical simulations of cluster lens- 



ing, it is currently not feasible to perform such simulations for suf- 
ficiently many combinations of the essential cosmological param- 
eters, i.e. the matter density parameter flo and the cosmological 
constant Therefore, conclusions can so far be drawn only for 
discrete points in the flo~flA plane. 

In a conceptually different approach, simple analytic, axially 
symmetric models have been used for describing the density pro- 
files of cluster lenses (see e.g. Wu & Mao 1996; Cooray, Quashnock 
& Miller 1999; Molikawa et al. 1999; Kaufmann & Straumann 
2000; Oguri, Taraya &. Suto 2001; Oguri et al. 2002; Molikawa 
& Hattori 2001). This method of investigation has the advantage 
that the computation of the probability for arcs satisfying a speci- 
fied property is fast and can easily be performed for a continuous 
and wide range of cosmological parameters, because the lensing 
properties of these models are perfectly known and fully described 
by analytic formulae. However, important effects like substructures 
or asymmetries in the matter distribution can at best be taken into 
account at an approximate level, and the correspondence between 
analytic and numerical models remains unclear. While the analytic 
studies by Cooray et al. (1999) and Kaufmann & Straumann (2000) 
find similar results as the numerical simulations regarding the sen- 
sitivity of arc statistics to the cosmic density parameter, their results 
are almost insensitive to the cosmological constant, in marked con- 
trast to Bartelmann et al. (1998), who found order-of-magnitude 
changes in the arc cross sections between low-density models with 
and without a cosmological constant. 

In this paper, we investigate whether the results of the analytic 
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and numerical approaches can be reconciled using a more realistic 
analytic lens model. Previous analytic studies of arc statistics com- 
monly used the singular isothermal sphere (hereafter SIS), which 
is computationally convenient, but has an unrealistic density pro- 
file and does not naturally reflect the theoretically expected and nu- 
merically demonstrated variations in halo concentration. Instead, 
we model cluster lenses as haloes with the density profile found by 
Navarro, Frenk & White (1997; hereafter NFW) in high-resolution 
simulations of haloes with a wide range of masses. The profile is 
flatter than isothermal near the core, and steeper outside, which 
leads to qualitatively and quantitatively different lensing proper- 
ties. We take the change of halo concentration with halo mass and 
cosmology into account, we distort the lensing potential elliptically 
in order to mimic cluster asymmetries, and we adapt the ellipticity 
by fitting numerical cluster models. We then compare the efficiency 
for producing arcs with a minimal length-to-width ratio of the ellip- 
tical NFW lenses with fully numerically simulated cluster lenses, 
and with the singular isothermal spheres for reference. 

The plan of the paper is as follows. In Sect. 2, we define the 
lensing cross section of both numerical and analytic models. In par- 
ticular, in Sect. 2.1 we describe the ray-tracing simulations, and 
we show the lensing cross sections produced by a sample of five 
galaxy clusters simulated in three different cosmological models. 
In Sect. 2.2, we deal with the lensing cross sections of axially sym- 
metric analytic models and discuss the differences between SIS and 
NFW models. In Sect. 2.3, we show how the lensing cross sec- 
tions change by assuming elliptical instead of circular sources. In 
Sect. 2.4, we discuss lensing by elliptically distorted NFW models 
and compare their arc cross sections with those of axially symmet- 
ric NFW lenses. In Sect. 3 we compare the analytic and the numer- 
ical models. Our conclusions are presented in Sect. 4. 



2 CROSS SECTIONS FOR LONG AND THIN ARCS 

The efficiency of a given lens for producing arcs with specified 
properties can be quantified by means of suitably defined lensing 
cross sections. They are defined as the areas on the source plane 
where sources must be located in order to be imaged as arcs satisfy- 
ing the required conditions on, e.g., size, shape, magnitude and the 
like. We focus in this paper on the cross sections of clusters for pro- 
ducing long and thin arcs, i.e. for arcs whose length-to-width ratio 
L/W exceeds a given minimum value {L/W)j^\j^. In the following 
subsections, we describe our method for computing strong lens- 
ing cross sections of nimierically and analytically modelled galaxy 
clusters. 

2.1 Numerical models 

For comparison to the analytic models, we compute the strong 
lensing cross sections of five numerically simulated cluster-sized 
dark-matter haloes, kindly made available by the GIF collabora- 
tion (Kauffmann et al. 1999). The same clusters have been used by 
Bartelmann et al. (1998). They were obtained from A''-body simu- 
lations performed in the framework of three different cosmological 
models. These are an Einstein-de Sitter model (hereafter SCDM); a 
flat, low-density universe with a matter density parameter = 0.3 
and a cosmological constant £1^ = 0.7 (ACDM); and an open, low- 
density model with f2() = 0.3 and f^A = (OCDM). The initial 
matter density in these models is perturbed about the mean accord- 
ing to a CDM power spectrum (Bond & Efstathiou 1984) with pri- 
mordial spectral index n = 1, normalised such that the local abun- 



dance of massive galaxy clusters is reproduced (e.g. Viana & Lid- 
die 1996). The complete list of cosmological parameters in these 
simulations is given in Bartelmann et al. (1998). The virial masses 
of the clusters at redshift z = range between ~ 5 x lO^^M© /h 
and~2x IO'^Mq/A. The masses of individual particles are 1.0 x 
IO^O/i-'Mq and 1.4 x 10^°h~^MQ for the high- and the Iow-Dq 
models, respectively. 

The lensing properties of these clusters are studied using a 
ray-shooting technique which was described in detail elsewhere 
(Bartelmann et al. 1998; Meneghetti et al. 2000, 2001). Here, we 
only briefly discuss some parameters used in these computations, 
and refer the reader to the cited papers for a more complete descrip- 
tion. 

Since strong lensing only occurs near the very cluster cen- 
tres, we select areas of \.5h~^ Mpc comoving side length centred 
on the lens centres. For all reasonable combinations of lens and 
source redshifts, this region is large enough to encompass the crit- 
ical curves of the cluster lens, close to which large arcs form. A 
bundle of 2048 x 2048 rays is shot through a regular grid on the 
lens plane, covering the studied region, and their paths to the source 
plane are traced. We consider different snapshots of the cluster sim- 
ulation at selected redshifts between zi = and zi = 1. For each 
snapshot, three projections of each simulated cluster on orthogonal 
planes are used for performing the lensing simulations. The typical 
angular resolution at z 0.3 is 0.18 arcsec, assuming the ACDM 
model. 

The sources are assumed to all he on one plane at redshift 
Zs = 1. Although real galaxy sources are obviously distributed in 
redshift, putting them all at the same redshift is acceptable because 
first, observations show that most of the sources which experience 
strong lensing effects are generally at redshifts near unity and sec- 
ond, the critical surface density changes very little with source red- 
shift, unless it is very close to the lens redshift. Moreover, we as- 
sume sources to be elliptical, with their axis ratios b/a uniformly 
distributed in the range 0.5-1.0. An initial set of sources is placed 
on a regular grid covering the source plane, and additional sources 
are added on sub-grids whose resolution is iteratively increased to- 
wards the lens caustics. Thus, sources are placed on an adaptive 
hierarchy of grids in order to improve the numerical efficiency of 
the method and to increase the probability for finding long and thin 
arcs. For correcting the statistics of the numerically simulated arcs, 
we assign to each source and all of its images a statistical weight 
w which is inversely proportional to the squared resolution of the 
sub-grid on which the corresponding source was placed. The finer 
the grid resolution, the lower is the statistical weight of its images; 
see also Eq. (1) below. 

Using the ray-tracing technique, we reconstruct the images of 
the background sources and measure their length, width and cur- 
vature radius. Our technique for image detection and classification 
was described in detail by Bartelmann & Weiss (1994) and adopted 
by Bartelmann et al. (1998) and Meneghetti et al. (2000, 2001). It 
results in a catalogue of simulated images which is subsequently 
analysed statistically. 

Each source is taken to represent a fraction of the source plane. 
The cells of the sub-grid with the highest resolution have area A, 
and the sources placed on its grid points are given a statistical 
weight of unity. The absolute lensing cross sections for a speci- 
fied image property are then determined by counting the statistical 
weights of the sources whose images have the required property. 
Specifically, we search for sources with a length-to-width ratio ex- 
ceeding a threshold (L/W)min- If a source has multiple images with 
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{L/W) > (L/W)nun, we multiply its statistical weight by the num- 
ber of such images. Therefore, the lensing cross section is 



--AY^WiWiUi 



(1) 



where Wt is unity if the ;-th source has images with (L/W) > 
(L/W)rom and zero otherwise, n,- is the number of images of the 
i-th source satisfying the required condition, and is the statisti- 
cal weight of the source. 



2.2 Axially symmetric models 

For sufficiently simple, axially symmetric lens models, the strong 
lensing cross sections can be also computed in an analytic or semi- 
analytic way. We consider two such models, the singular isothermal 
sphere (SIS) for reference, and haloes with the NFW density profile 
(Navarro et al. 1997). The SIS has been very widely used in many 
previous studies on arc statistics because of its computational sim- 
plicity. The NFW density profile fits the results of highly resolved 
numerical halo simulations and thus provides a much more realistic 
description of cluster density profiles than the SIS profile. We focus 
on the NFW profile here, but give some results for SIS lenses for 
later comparison. 

The SIS density profile is given by 



p(r): 



(2nG;-2) ' 



(2) 



where is the velocity dispersion. The NFW density profile is 
Ps 



P(^) 



(r/rs)(H-r/rs)2 



(3) 



where ps and are characteristic density and distance scales, re- 
spectively (see Navarro et al. 1997). These two parameters are not 
independent, but related to a single parameter, which can be taken 
as the cluster mass. It is important to note that Ps and also depend 
on the cosmological model, hence the lensing properties of haloes 
with identical mass are different in different cosmological mod- 
els, a crucial property which the simpler singular isothermal sphere 
models do not share. The NFW halo falls off more steeply than 
isothermal at radii beyond r^, but flattens towards the halo centre. 
These different features lead to markedly different lensing prop- 
erties of the NFW compared to the SIS model, as will be shown 
below (see also the discussion in Perrotta et al. 2002). 

For axially symmetric lens models, the lens equation is essen- 
tially one-dimensional. Define the optical axis as the line running 
from the observer through the lens centre, and introduce physical 
distances from the optical axis ^ and r\, respectively, on the lens 
and source planes. Fixing a length scale on the lens plane, we 
can define the dimensionless distance x = ^/^o- The chosen length 
scale ^0 is projected onto the length scale r|o = i,oDs/D\ on the 
source plane, where Di and are the angular-diameter distances 
to the lens and source planes, respectively. In analogy to the lens 
plane, a dimensionless distance from the optical axis y = r\/r\Q can 
now be defined on the source plane. The lens equation, relating the 
position of an image on the lens plane to that of its source on the 
source plane is then 



y-- 



-alx) 



(4) 



where a{x) is the reduced deflection angle at distance x from the 
lens centre, caused by the lensing mass distribution. It is the gradi- 
ent of the lensing potential \|/, 



a{x) ■■ 



di|/(x) 
dx 



(5) 



See, e.g. Schneider, Ehlers & Falco (1992) and Narayan & Bartel- 
mann (1997) for more details. 

Local imaging properties are described by the Jacobian ma- 
trix of the lens mapping. It has two eigenvalues. At and A.t, which 
describe the image distortion in the radial and the tangential direc- 
tions, respectively. For an axially symmetric lens, they are 

ArW = l-^, AtW = l-"- (6) 

Analytic expressions for the lensing potential at any dis- 
tance x are easily obtained for the SIS and the NFW lens models. 
Therefore, the eigenvalues of the Jacobian matrix can straightfor- 
wardly be computed for both lens models considered here. 

For a SIS lens, the lensing potential is 

¥W = \x\ , (7) 

if ^0 = 47t(cyi,/c)^DiDis/Ds is chosen as a length scale, where Di^ 
is the angular diameter distance between lens and source. Hence, 
the deflection angle and the eigenvalues are 

1 



a{x) -■ 



X 

M 



1 



At = 1- 



(8) 



For an NFW lens, we take = ''s for simplicity and deflne 
Ks = psrsE~\ where £cr = [c^/(47iG)] [Ds/{DiDis)] is the critical 
surface mass density for strong lensing. We then have the lensing 
potential 



V|/(x) = 4Ks 



- -2arctanh^ 
2 2 



1+x 



which impUes the deflection angle 



a{x) -■ 



4Ks 

X 



2 , /l-^ 

In — \ , arctanh^ / 

1+x 



(9) 



(10) 



from which the eigenvalues A^-^t can straightforwardly be derived. 
Several different aspects of lensing by haloes with NFW or gener- 
aUsed NFW profiles can be found in Bartelmann (1996), Wright & 
Brainerd (2000), Li & Ostriker (2002), Wyithe, Turner & Spergel 
(2001), Perrotta et al. (2002). It can easily be verified that the po- 
tential [Eq. (9)] satisfies the Poisson equation V^v|/ = 2k, with k 
given in Bartelmann (1996). 

It is an important feature of the NFW lensing potential 
[Eq. (9)] that its radial profile is considerably less curved near the 
centre than the SIS potential [Eq. (7)]. Since the local imaging 
properties are determined by the curvature of \\f, this immediately 
implies substantial changes to the lensing properties. At fixed halo 
mass, the critical curves of an NFW lens are closer to its centre 
than for a SIS lens because of its flatter density profile. There, the 
potential is less curved, thus the image magnification is larger and 
decreases more slowly away from the critical curves. Therefore, 
NFW lenses are less efficient in image splitting than SIS lenses, but 
comparably efficient in image magnification. What is more impor- 
tant here is that any additional shear added to a flat lensing poten- 
tial (e.g. by asymmetries) much more strongly extends the critical 
curves than for a steeper potential. We will see the consequences 
below. 

The inverse of the eigenvalues Ar.t of the Jacobian matrix gives 
the radial and the tangential magnifications at the radial distance x 
from the lens centre. The points satisfying the conditions Ar,t = 
form the radial or tangential critical curves, respectively, where the 
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corresponding magnifications tend to infinity. For SIS lenses, 7^ is 
always unity, hence the radial critical curve does not exist in this 
case, so the images are not radially magnified. The tangential crit- 
ical curve is the circle |x| = 1. For NFW lenses, the critical curves 
have to be found numerically. In the limit of small Ks, when the 
critical curves are close to the lens centre, we find to first order in x 



^ 2exp 



1+Ks 

2Ks 



At 
2.718 



(11) 



By means of the lens equation, the critical curves are mapped 
to the corresponding caustic curves in the source plane. For the 
axially symmetric lenses we are considering, the tangential caustic 
degenerates to the point where the optical axis intercepts the source 
plane. 



2.3 Elliptical sources 

The length-to-width ratio of an image depends on the local lens 

properties and on the ellipticity of the sources. We derive and com- 
pare results for both circular and elliptical sources. Since the source 
galaxies are typically much smaller than the length scale on which 
the lens properties change substantially, we can assume that the 
length-to-width ratio of the images of a circular source is given by 
the ratio of the two eigenvalues q\ = 'k^ /Xi at the image position. 

An elegant model for introducing elliptical sources into ana- 
lytic lensing calculations was recently proposed by Keeton (2001). 
Following his approach, the observed length-to-width ratio LjW 
is a function of three variables, which are the intrinsic axis ratio 
= a/b and the intrinsic position angle 6 of the source, and the 
eigenvalue ratio q\ of the lens at the image position: 

1/2 

(12) 



L/W = 

where 
T = 

Q = 



T + (T^ 



r-(r2-4e)i/2 



) cos^ e , 



(13) 
(14) 



The length- to- width ratio L/W of each image can be related 
to the source position y by means of the lens equation. For sources 
with fixed intrinsic shape and orientation lensed by axially symmet- 
ric lenses, L/W monotonically decreases as the source is moved 
from y = (where the degenerate tangential caustic is located) 
away from the lens centre. Thus the lensing cross section for arcs 
with L/W > {L/W)jtim from sources with intrinsic axis ratio and 
position angle 6 can be written as 



'(i/W)„ 



„(9s,e): 



(15) 



where y^in is the distance from the optical axis where L/W falls 
below (L/W)^i„. 

In order to compute the final lensing cross section, we have to 
average these cross sections over all the possible values of and 6- 
For doing so, we assume that 9 and q^^ are uniformly distributed 
in the ranges [0,7l] and [0.5, 1], respectively. This agrees with the 
assumptions underlying the numerical simulations. Therefore, the 
final lensing cross section is given by 



-/7 

K Jl Jo 



0(L/w)^(?s,e)de^^ 



(16) 



Both SIS and NFW lenses can produce multiple images for 
suitable source positions. For SIS lenses, multiple images occur 
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Figure 1. Cross sections 07.5 for the formation of arcs with length-to-width 
ratio larger than 7.5 as functions of lens redshift. Results are shown for a 
virial lens mass of M = lO'^ in a ACDM cosmology. Two different 

lens models are used: the axially symmetric NFW lens model (solid curves), 
and the singular isothermal sphere for comparison (dashed lines). Thick 
and thin curves were obtained assuming intrinsically circular and elliptical 
sources, respectively. 



only for sources lying inside the Einstein ring, which has the an- 
gular radius 9e = {^na^/c^) (Di^/Ds). If this is the case, the lens 
produces an arc and a counter-arc. As the source moves towards 
the Einstein ring, the tangential magnification of the counter-arc 
decreases and tends to zero. Then, the counter-arc disappears (see 
Narayan & Bartelmann 1997 for further detail). The main "arc" re- 
mains, but its tangential magnification goes to unity as the source 
moves away from the optical axis towards infinity. 

For NFW lenses, multiple images occur only for sources with 
y < yc,T, where yc,T is the radius of the radial caustic. As shown by 
Bartelmann (1996), the NFW density profile always has a radial 
critical curve (and thus a radial caustic) for any combination of Ps 
and rs. Sources inside the radial caustic therefore always produce a 
radial image, an arc, and a counter-arc. We neglect the radial image 
because it is radially oriented. As for SIS lenses, the arc and the 
counter-arc have different length-to- width ratios L/W. In fact, the 
counter-arc is radially magnified as the source moves towards the 
radial caustic, and its length-to-width ratio decreases more rapidly 
than that of the main tangential arc. Therefore, two separate cross 
sections must be computed for both SIS and NFW lens models, one 
for the arc and one for the counter-arc, and the total lensing cross 
section is the sum of these two contributions. 

To give an example. Fig. 1 shows the lensing cross section for 
arcs with length-to-width ratio L/W > 7.5 produced from sources 
at redshift Zs = 1 by a lens with mass W^^h-^MQ at redshifts 
between zero and unity. Results are given for the ACDM model. 
The solid lines correspond to NFW lenses, the dashed lines to SIS 
lenses. Thick and thin lines indicate the results for elliptical and 
circular sources, respectively. We use here the length-to-width ra- 
tio threshold of 7.5 instead of 10 in order to improve the statistics 
of numerically simulated arcs, and to reduce corresponding fluctu- 
ations in the cross sections. 

As expected, the lensing cross section of the NFW model is 
smaller than for SIS lenses because of its flatter density profile. In 
particular, NFW lenses lose their strong-lensing efficiency when 
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the lens approaches the observer or the sources. Conversely, due 
to its unrealistically steep and scale-free density profile, the SIS 
remains an efficient strong lens even when it is located very close 
to the observer or to the source. 

The relative increase of the cross section for the intrinsically 
elliptical sources is virtually independent of the lens redshift, but it 
does depend on the lens model. In fact, the relative change in the 
lensing cross sections for NFW lenses is approximately twice as 
high as for SIS lenses, quite independent of the lens redshift. This 
is again due to the shallower lensing potential of NFW compared 
to SIS lenses. 



2.4 Elliptical models 

The construction of lens models with elliptical or pseudo-elliptical 
isodensity contours is generally quite complicated (see e.g. Kas- 
siola & Kovner 1993; Kormann, Schneider & Bartelmann 1994; 
Golse & Kneib 2002). Starting from an elliptical lensing potential 
is computationally much more tractable, but has the disadvantage 
that the mass distribution corresponding to the elliptical potential 
can become dumbbell-shaped even for moderate ellipticities, which 
is unwanted for galaxy lenses. Galaxy clusters, however, are less re- 
laxed and exhibit substructure, so for them dumbbell-shaped mass 
distributions are uncritical. For this reason, we prefer constructing 
elliptical cluster lens models starting directly from the effective 
lensing potential. Moreover, since the NFW density profile gives 
a much more realistic reproduction of the simulated clusters than 
the SIS, we chose to generalise only the NFW lens model to the 
elliptical case. 

The lensing potential of an axially symmetric NFW lens was 
given in Eq. (9) above. We introduce the ellipticity e= \ —b/a, 
where a and b are the major and minor axes of the ellipse, by sub- 
stituting 



«2 



(1- 



+ x^(I-e) 



(17) 



where xi and X2 are the two Cartesian components of x, x^ = Xj + 
x^ - This ensures that the mass inside circles of fixed radius remains 
constant as the ellipticity changes. 

The Cartesian components of the deflection angle are then 

dvf _ Xi 
3x1 



3\|/ 

3x2 



X2(I -e) 



a{x) 



(18) 



a{x) 



where a(x ) is the unperturbed (i.e. axially-symmetric) deflection 
angle at the distance X from the lens centre. 

Using these formulae, deflection-angle fields for different val- 
ues of the ellipticity e are readily computed. Some examples for 
deflection-angle maps are displayed in Fig. 2. Obviously, the shape 
of the contours becomes more and more elliptical as e is increased. 
We analyse the lensing properties of these deflection-angle fields 
using the ray-tracing technique and, finally, compute the lensing 
cross sections, adopting the same techniques as applied to the nu- 
merical models. 

Increasing the ellipticity of the lensing potential strengthens 
the shear field of the lens, and consequently the tangential caustic 
expands and changes. Examples for the change of the caustics with 
ellipticity e are shown in Fig. 3, which refers to a halo with mass 
M = lO'^ h^^ Mr.) at redshift z = 0.3, and the underlying cosmol- 
ogy is the ACDM model. As discussed above, the radial and tan- 
gential caustics are a circle and a point, respectively, for the axially 
symmetric models. Increasing e, the caustics stretch, develop cusps, 
and enclose a growing area. Thus the lensing cross sections are ex- 
pected to grow rapidly. This is confirmed by Fig. 4, which shows 
cross sections for arcs withL/lV > 7.5. These results were obtained 
through ray-tracing simulations in which deflection angle maps for 
two haloes of mass M = 7.5 x lO''* h'^ Mr, and M = lO'^ A"' Mq 
were used. Again, we adopt the ACDM model universe and put 
the lens at redshift z = 0.3. Increasing the ellipticity of the lensing 
potential from e = to e = 0.5, the cross section increases approx- 
imately by a factor of 30, almost independently of the lens mass 
considered. 



3 COMPARISON OF ANALYTIC AND NUMERICAL 
MODELS 

We can now compare the strong lensing cross sections of the nu- 
merical and analytic models introduced in the preceding sections. 
For that purpose, we focus on the lensing cross sections for the for- 
mation of arcs with length-to-width ratios L/W larger than 7.5 and 
10. As described before, we compare them to the lensing properties 
of five numerical models of galaxy cluster haloes, picking simula- 
tion snapshots at twelve different redshifts between zero and unity. 
For clarity, we present here the results obtained for the most mas- 
sive halo only. The behaviour of the cross sections for the other 
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Figure 3. Caustic curves produced by an NFW halo with mass M = 10'^ h- ' Mrj at redshift z = 0.3 in a ACDM universe. Different panels show results for 
different ellipticities e of the lensing potential: e = 0, i.e. axially symmetric (left panel); e = 0.2 (central panel) and e = 0.4 (right panel). The figure demonstrates 
that even small or moderate elliptical distortions substantially stretch the caustic curves, and thus cause the arc cross sections to grow considerably. 



numerical models is in good qualitative and quantitative agreement 
with that obtained for this cluster. 

The results are illustrated in Fig. 5, where the dotted lines refer 
to the fully numerically simulated cluster, while solid and dashed 
lines represent the cross sections of NFW and SIS lenses, respec- 
tively, having the same virial mass as the numerical cluster model. 
Finally, the shaded regions in the same plots show the cross sections 
obtained by elliptically distorting the NFW lensing potential with 
ellipticities e in the range between e = 0.2 and e = 0.4 (lower and 
upper limits, respectively). Results are shown for L/W = 7.5 (up- 
per panels) and L/W = 10 (lower panels), and they were obtained 
for the SCDM (left panels), ACDM (central panels), and OCDM 
models (bottom panel). 

First, we checked the calibration of the analytic relative to the 
numerical cross sections. For doing so, we performed ray-tracing 
simulations using the axially symmetric deflection angle maps for 
the NFW lens model, and used them for determining the strong- 
lensing cross sections. The results are shown as filled dots in the 
panels of Fig. 5. The very good agreement with the analytic esti- 
mates obtained as explained in Sect. 2.2 demonstrates the reliability 
of our numerical technique. 

The general trends in the lensing cross sections shown in Fig. 5 
can be understood as follows. The strong-lensing efficiency of a 
mass distribution depends on several factors. First, for the light 
coming from the sources to be focused on the observer, the lens 
must be located at a suitable distance from both the observer and 
the sources. Second, the larger the (virial) mass of the lens is, the 
stronger are the lensing effects it produces close to its centre. Fi- 
nally, the more concentrated the lens is, the thinner are the long 
arcs expected to be. 

Since the lens mass grows with decreasing redshift because 
new material is accreted by the halo and deepens its potential well, 
the lensing cross sections are expected to grow as well. On the other 
hand, when the lens is too close to the observer, the cross section is 
geometrically suppressed, unless the lens surface density profile is 
sufficiently steep and scale-free, as for the SIS lenses. In fact, in this 
case the focusing by the lens is strong enough to allow observers 
to see strongly distorted images of background sources also in very 
near lenses, i.e. at relatively small redshifts (see also Fig. 1). 

Moreover, the results are expected to depend on cosmology. In 
fact, in the OCDM and ACDM models, galaxy clusters become ef- 
ficient lenses at higher redshift than in the SCDM scenario, because 
they form earlier and are centrally more concentrated. This is con- 




0,1 0,2 0.3 0,4 0,5 
ellipticity 



Figure 4. Strong-lensing cross section 07 5 of elliptically distorted NFW 
lenses for arcs with length-to-width ratio larger than 7.5 as a function of the 
ellipticity e of the lensing potential. Tick and thin lines were calculated for 
haloes ofmassM= lO'^ft^'M., andM = 7.5 x 10'"*/!"' Mq, respectively, 
both placed at z = 0.3 in a ACDM cosmology. 



firmed by our numerical results. In the Einstein-de Sitter case both 
CT7.5 and <Jio are completely negligible for redshifts z 0.4-0.5 and 
peak at z ~ 0.2, while for low-density models the cross sections are 
largest in the redshift range 0.3 ^ 0.5. A more complete dis- 
cussion on the dependence of lensing cross sections on cosmology 
can be found in Bartelmann et al. (1998). 

In this paper, we concentrate on comparing the cross sections 
of the analytic and numerical models. As Fig. 5 shows, the nu- 
merical models generally have much larger cross sections than the 
analytic models. In particular, the cross sections for axially sym- 
metric NFW lenses are almost two orders of magnitude smaller, 
quite independent of the cosmological models considered. For SIS 
lenses, the discrepancy with the numerical models is only partially 
compensated by the unrealistically steep central density profile, but 
the estimated values of 07 5 and <Jio remain too low. Introducing 
the elliptical distortion into the NFW lens model allows the cross 
section to increase by roughly an order of magnitude compared to 
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Figure 5. Lensing cross sections for arcs witli lengtli-to-width ratio larger tlian 7.5 (07 5; top panels) and 10 (Cio; bottom panels) as a function of the lens 
redshift. Results for different cosmologies are shown: SCDM (left panels), ACDM (central panels) and OCDM (right panels). Dotted lines show results 
obtained from the analysis of the most massive halo in the numerically simulated sample; solid and dashed lines represent the cross sections predicted by the 
axially symmetric NFW and SIS models, respectively, with the same viiial mass as the numerically modelled halo. The shaded regions mark the cross sections 
obtained with the elHptically distorted NFW lens model, with ellipticities ranging between e = 0.2 and e = 0.4 (lower and upper Hmits, respectively). The 
filled dots are the lensing cross sections found by ray-tracing simulations with the deflection angle maps of the axially symmetric NFW lens model. Their 
agreement with the analytic curves demonstrates the reliability of our numerical method. 



the axially symmetric NFW model, but even then the analytic cross 
sections fail to reproduce the numerical cross sections unless unre- 
alistically high values of e are adopted. 

Fig. 5 also demonstrates that apart from the cross-section am- 
plitude, the analytic models miss another important feature of the 
numerical results, which show steep increases and decreases re- 
flecting merger events. While a merging sublump approaches a 
cluster, the cross section tends to increase because of the increasing 
tidal field, and as the cluster relaxes following the merger event, the 
cross section decreases again. It becomes quite clear from Fig. 5 
that such events play an important role in understanding realistic 
cluster cross sections, and they cannot reasonably be captured in 
analytic models unless their prime advantage of being computa- 
tionally fast is sacrificed. 

3.1 EUipticity estimates 

It is important to check which ellipticities are typical for the lensing 
potential of the numerical clusters. To this end, we compare the 
deflection-angle maps constructed for the numerical haloes with 
those obtained for elliptical NFW lenses with identical virial mass. 
We estimate the best-fit ellipticity e by minimising the mean-square 
deviation 



where a,- and are the deflection angles of the light ray passing 
through the ;-th grid point on the numerical and the elliptical lens 
model, respectively. The summation is done over all grid points 
contained in the central region of roughly ~ 200h^^ kpc comoving 
side length. This way, we measure the ellipticity of the lensing po- 
tential in the central region close to the critical curves which is the 
most relevant for our purposes. For this test, we use the simulation 
snapshots at z = 0.27, which is approximately the redshift where 
the numerical haloes reach their maximum lensing efficiency. The 
complete sample of five numerical clusters and three independent 
projections per cluster was used for this analysis. 

The results are reported in Tab. 1. Median ellipticities of 
e ~ 0.3 are found, with only a weak dependence on cosmology. 
Indeed, only little evolution of the ellipticity is found going from 
low to high-density cosmological models. The median ellipticities 
in the OCDM and the ACDM models are slightly smaller than in 
the SCDM case, as expected, given that clusters form earlier in 
these models and have more time to relax. Anyway, as the semi- 
interquartile ranges (SIQR) show, the ellipticity distributions are 
quite broad, in particular for SCDM. The same analysis has also 
been repeated using other simulation snapshots, corresponding to 
redshifts in the range 0.2 ^ z ^ 0.7. The results are quite similar 
to the previous ones and therefore not given in the table. The typi- 
cal median of e is around 0.3, and exceeds 0.4 in only a few cases. 
In conclusion, the ellipticities measured in the lensing potential of 
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Table 1. EUipticities e of the lensing potential obtained by fitting the 
deflection-angle maps of the eUiptically distorted NFW lens model to fully 
numerical deflection angles. We give the medians (second column) and 
semi-interquartiles ranges (SIQR; third column) for the ellipticities found 

in our cluster sample for different cosmological models (indicated in the 
first column). Results are shown for the simulation snapshots at redshift 
z = 0.27. 



cosmological model 


median 


SIQR 


SCDM 


0.320 


0.102 


ACDM 


0.305 


0.032 


OCDM 


0.270 


0.052 



the numerical clusters are insufficient for reconciling the cross sec- 
tions of eUiptically distorted NFW models with the fully nimierical 
results. 

3.2 Substructure estimates 

The remaning difference between numerical and analytical cross 
sections must be attributed to some factors which are missing from 
the analytical models. The most important of those is certainly the 
presence of substructure in the lensing mass distribution. 

Deviations of the projected mass distribution of a numerically 
simulated cluster from the predictions of circular or elliptical mod- 
els can be quantified by means of a multipole expansion of its sur- 
face density field. For performing this analysis, we first define a ref- 
erence frame whose origin coincides with the surface density peak 
of the cluster. Then, starting from the particle positions in the N- 
body simulations, we compute the surface density at discrete radii 
r„ and position angles ^i^, with r„ and (|)^ taken from the intervals 
[0, 1.5]/!^'Mpc and [0,27t], respectively. For any given r„, each dis- 
crete sample of data £(r„ , ^i^) can be expanded into a Fourier series 
in the position angle, 

Arn,ifk) = E 5,(r„)e-'* , (20) 
where the coefficients Si{r„) are given by 

Si{r„) = £ E(r„,0^)e'* , (21) 

and can be computed using fast-Fourier techniques. We define the 
power spectrum /*„(/) of the multipole expansion / as /"«(/) = 

Axially symmetric and elliptical models have very simple 
multipole expansions. For a circular mass distribution, only the 
monopole, 1 = 0, contributes to the sum in Eq. (20). An elliptical 
mass distribution has one more contribution from the quadrupole, 
1 = 2. The dipole term, Z = 1, is zero for both these models. Of 
course, the surface density fields of numerically simulated clusters 
are much more complex than those of our analytic models. Their 
multipole expansions contain a dipole and also multipoles of higher 
order than two, which correspond to substructures spanning open- 
ing angles of order ~ n/l when seen from the cluster centre. 

Therefore, in order to quantify the amount of substructure and 
the degree of asymmetry in the mass distributions of our numeri- 
cally simulated lenses at any distance r„ from the cluster centre, we 
can use the power spectra P„{t). In particular, we define an inte- 
grated power Pint (?■«). which is the sum of the power spectral den- 
sities of the dipole and of all multipoles of higher order than two, 



i.e. we subtract the monopole and quadrupole contributions from 
the total integrated power, 

^'int('-«) = £ Pn{l)-Pn{0)-Pn{2) . (22) 

This quantity measures the deviation from an elliptical distribution 
of the surface mass density at a given distance r„ from the clus- 
ter centre. In order to suppress the dependency on the radial coor- 
dinate, we also compute the area-weighted averaged value of Pint 
inside radius r as 

^int('-) = 4 E Pi^t{rn)r„Ar„ . (23) 

rn<r 

In Fig. 6, we show how the integrated power Pint('') changes as 
a function of the cluster redshift for those lenses whose cross sec- 
tions were plotted in Fig. 5. Again, we show the results for all three 
cosmological models. For making the integrated power at differ- 
ent radii comparable, we normalize them to the power of the cor- 
responding monopole. Moreover, we average the results obtained 
for the three projections of the same cluster. The heavy dotted and 
dashed lines show the results for radii ri = 1.25/!~^Mpc, compa- 
rable to the cluster virial radius, and = 155/i^'kpc. For com- 
parison, the heavy solid curves indicate the area-weighted power 
inside circles of radius n as a function of redshift. The small val- 
ues reached by the heavy dashed curves indicate that in the very 
central region of the cluster the mass distribution is dominated by 
the monopole and the quadrupole terms, i.e. in the innermost re- 
gions of the clusters, the surface mass density has elliptical iso- 
contours. Moreover, these curves are flat, which means that the sur- 
face density contours remain elliptical at all redshifts between zero 
and unity. On the other hand, as shown by the heavy dotted and 
solid lines, at distances comparable to the virial radius or smaller 
the contribution to the power from the dipole and the multipoles 
of higher order than two is large and can even exceed 15% of the 
monopole contribution. 

A quick comparison of these curves to the lensing cross sec- 
tions in Fig. 5 shows that the redshifts where the contributions of 
the dipole and the higher-order multipoles are largest correspond 
quite well to those where the numerical cross sections deviate most 
strongly from those of the elliptical models. 

This is most obvious between redshifts 0.2 and 0.4, where also 
the geometrical lensing efficiency is largest. At lower and higher 
redshifts, the lensing cross sections drop for two major reasons; 
first, the lenses are too close to the observer or to the sources; 
second, the cluster virial mass decreases with increasing redshift. 
For better comparing the dependences on redshift of the integrated 
power and the lensing cross sections, we rescale Pint with the ef- 
fective lensing distance, D = D\D^slDf, and with the virial cluster 
mass. The thin lines in Fig. 6 show Pjnt after rescaling. The damp- 
ing effect almost completely removes the peaks in the integrated 
power at redshifts higher than z ~ 0.5 — 0.6, as well as those at red- 
shifts lower than z ~ 0.2, and make the curves much more similar 
to the curves displaying the lensing cross sections as a function of 
redshift in Fig. 5. 

This correlation between the higher-order multipoles of the 
mass distributions and the cross sections proves that substructures 
and aymmetries are the dominant reason for the discrepancy be- 
tween numerical and analytical lensing cross sections. We verified 
that the peaks in the integrated power Pint occur when clumps of 
matter enter within the virial radius of the respective clusters. Then, 
the lensing cross section grows for two reasons; first, because the 



© 0000 RAS, MNRAS 000, 000-000 



Cluster cross sections for strong lensing: analytic and numerical lens models 9 






Figure 6. Integrated multipole power as a function of redshift for the most massive cluster in our numerically simulated sample. Results for three different 
cosmologies are shown; SCDM (left panel). ACDM (central panel) and OCDM (right panel). Heavy dotted and dashed lines show the results of multipole 
expansions of the cluster surface mass density at cluster-centric distances of ri = 1.25/i^'Mpc and rg = 155/!"'kpc, respectively, and heavy solid lines show 
the area-weighted averages of the integrated power within ri (see text for more details). Light lines show the results of rescaling the heavy curves with the 
cluster virial masses and the effective lensing distance. 



cluster mass increases, and second, because of the increased shear 
produced by the substructures falling towards the cluster centres. 



4 SUMMARY AND DISCUSSION 

In this paper, we have introduced an analytic lens model for de- 
scribing strong lensing by galaxy clusters. It improves upon pre- 
viously used models in four ways. First, the commonly adopted 
singular isothermal density profile is replaced by the NFW profile, 
which is much more adequate for cluster-sized haloes. Second, this 
model allows the effect of cosmology on the lens concentration to 
be included. Third, we elliptically perturb the model for approxi- 
mating the substantial effect of cluster asymmetries on their strong- 
lensing cross sections. Fourth, we take into account that sources 
are intrinsically elliptical rather than circular. We adapt the ellip- 
ticity of this lens model to fit numerically simulated clusters and 
compare the strong-lensing cross sections of the analytic model to 
those of fully numerically simulated clusters. Results for singular 
isothermal sphere lenses are given for comparison. Our results can 
be summarised as follows: 

• The cross section of the axially symmetric NFW lens model 
for arcs with length- to- width ratio larger than 7.5 and 10 are almost 
two orders of magnitude smaller than those of simulated clusters. 

• The axially symmetric singular isothermal lens model is more 
efficient for strong lensing because its density profile is steeper in 
the core, but its large-arc cross sections are still smaller than the 
numerical ones by typically an order of magnitude, except at low 
redshift where both the numerical models and the NFW model fail 
in producing any lensing effect. 

• The comparatively flat lensing potential of haloes with the 
NFW density profile makes NFW lenses poor image splitters, but 
efficient magnifiers. In addition, the flat potential renders the lens- 
ing properties very sensitive to changes in the tidal field, because 
small deformations of the potential can lead to large shifts of the 
critical curves. 

• The cross sections of our new elliptical model increase steeply 
and monotonically as a function of ellipticity. An ellipticity of 
~ 0.3 typically increases the cross sections by an order of mag- 
nitude. However, for the model to reproduce the arc cross sections 



of the numerical clusters, unrealistically high ellipticities e <; 0.5 
are required. 

• Comparing the deflection angle maps to those of elliptical 
NFW lens models with variable e, we estimated the ellipticity of 
the lensing potential in the central region of the numerically mod- 
elled haloes and found typical ellipticities of e --^ 0.3, substantially 
below values required to solve the discrepancy between analytic 
and numerical lens models. 

• The change of the fully numerical cluster cross sections with 
time exhibits pronounced signatures of merger events. As merging 
sublumps approach the cluster centre, their tidal field markedly in- 
creases the strong-lensing cross section. We have verified this by 
means of a multipole expansion of the cluster surface density field, 
which shows that larger deviations of the lensing cross section of 
the numerical clusters from the prediction of the analytic models 
arise when the contribution to the surface density power spectra 
from the dipole and from higher-order multipoles of high order. 
Our comparison of elliptical analytic lens models with fully nu- 
merical models shows that an adequate description of such events 
is necessary for an accurate calculation of arc cross sections. 

We conclude that even our improved analytic model is unable 
to reproduce the strong lensing properties of realistic cluster mod- 
els, which we assume the fully numerically simulated haloes to be. 
We have seen that the axially symmetric NFW lens model under- 
predicts the number of arcs with length-to-width ratio exceeding a 
given threshold by approximately two orders of magnitude com- 
pared to the fully numerical results. 

At present, the NFW density profile can be considered the 
most realistic model profile for cluster haloes (see, however, Moore 
et al. 1999; Jing & Suto 2000). Moreover, it allows to take the effect 
of varying halo concentration on the strong lensing efficiency into 
account, which is a substantial advantage compared to the com- 
monly used SIS density profile. Finally, the profile permits analyt- 
ical calculations of the relevant lensing properties. Therefore, the 
NFW density profile appears ideally suited for constructing ana- 
lytic models for strong lensing by clusters, in particular if elliptical 
distortions are included, as we have done in our extension of the 
model. 

However, the differences remaining between the fully numer- 
ical and analytic approaches indicate that analytic modelling is still 
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insufficient for properly and accurately describing strong lensing Wu, X.-R, Mao, S., 1996, ApJ, 463, 404 
by galaxy clusters. We showed that the most important missing WyitheJ.S.B.,TumerE.L.,SpergelD.N., 2001, ApJ, 555, 504 
factors are the presence of substructures within the clusters, and 
the tidal field of the surrounding matter distribution. In fact, signif- 
icant substructure is abundant in and around numerically simulated 
haloes. They enhance the shear field around the clusters, increas- 
ing the length of the critical curves and consequently increase the 
probabiUty of forming long arcs. 

Of course, the elliptically distorted NFW lens model is ade- 
quate for good qualitative calculations of arc probability, as Fig. 5 
shows. However, in order to derive precise constraints on cosmol- 
ogy or cluster structure and evolution, more realistic cluster lens 
models are required, for which numerical simulations seem to be 
the only reliable choice. 
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